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ABSTRACT 

We report a statistical study of flares observed with the Soft X-ray Telescope 
(SXT) onboard Yohkoh in the year of 2000. We measure physical parameters 
of 77 flares, such as the temporal scale, the size, and the magnetic flux density 
and find that the sizes of flares tend to be distributed more broadly as the GOES 
class becomes weaker and that there is a lower limit of magnetic flux density that 
depends on the GOES class. We also examine the relationship between these pa- 
rameters and find weak correlation between temporal and spatial scales of flares. 
We estimate reconnection inflow velocity, coronal Alfven velocity, and reconnec- 
tion rate using above observed values. The inflow velocities are distributed from 
a few km s _1 to several tens km s _1 and the Alfven velocities in the corona are 
in the range from 10 3 to 10 4 km s _1 . Hence the reconnection rate is 10~ 3 - 10~ 2 . 
We find that the reconnection rate in a flare tends to decrease as the GOES 
class of the flare increases. This value is within one order of magnitude from 
the theoretical maximum value predicted by the Petschek model, although the 
dependence of the reconnection rate on the magnetic Reynolds number tends to 
be stronger than that in the Petschek model. 

Subject headings: Sun: corona — Sun: flares — Sun: magnetic fields 

1. INTRODUCTION 

In the energy release process in solar flares, magnetic reconnection is generally considered 
to play a key role. The reconnection rate is an important quantity, because it puts critical 
restrictions on the reconnection model. It is defined as Ma = Via/V\ in non-dimensional 
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form, where V- ni is the velocity of the reconnection inflow and Va is the Alfven velocity. It 
gives the normalized value of the reconnected flux per unit time. In spite of its importance, 
what determines the reconnection rate in flares is still a question. In the steady reconnection 
model of Sweet (1958) and Parker (1957), the reconnection rate is Ma = (-R m ) -1 / 2 , where 
Rm = (VaL/t]) is the magnetic Reynolds number defined with the Alfven velocity and rj is 
the magnetic diffusivity. In the solar corona, if the resistivity is attributed to the Coulomb 
collision (Spitzer 1956), the typical value of R m is R m ~ 10 14 , which means M A ~ 1CT 7 and 
estimated timescale of flare is about one year. It is, of course, too slow to explain flares whose 
timescale is about 10 2 — 10 3 seconds. On the other hand, Petschek (1964) pointed out that the 
previous model lacks the effect of wave and suggested his model with M A < 7r/(8m(8i? m )) 
. The special feature of this model is that Ma has weak dependence on R m . In this model 
Ma < 1CT 2 when R m ~ 10 14 and the estimated timescale is consistent with observed value 
(see also Nitta 2004). 

The Sweet-Parker model and the Petschek model shown above are steady reconnection 
models. The physics of nonsteady reconnection is not established clearly. Shibata & Tanuma 
(2001) discussed the plasmoid-induced-reconnection model. Since flares are nonsteady phe- 
nomena, they suggested possibility that the reconnection rate in a flare is determined by the 
nonsteady and dynamic processes, such as plasmoid ejections. Another uncertainty in the 
model of the reconnection mechanism is uncertainty of the magnetic Reynolds number R m . 
It was suggested that finite levels of MHD turbulence can lead to increase of the magnetic 
diffusivity, i.e., decrease of the effective magnetic Reynolds number (Matthaeus & Lamkin 
1985; Tajima & Shibata 1997, §3.3). 

Though there are some observational case studies about the reconnection rate, statis- 
tical studies are very few so far. Dere (1996) examined the spatial and temporal scales of 
reconnection events in the solar atmosphere described in the literature and found that the 
reconnection rate is on the order of 0.001 - 0.1. Tsuneta (1996) analyzed the temperature 
structure of 1992 February 21 LDE flare in detail and found that the reconnection rate is 
0.07. He assumed isothermal switch-off shock and used the condition of the shock to derive 
the reconnection rate. The flare on 1992 January 13 was analyzed by Tsuneta et al. (1997) 
and the reconnection rate in this flare was estimated at 0.06. Ohyama & Shibata (1997) 
and Ohyama & Shibata (1998) studied 1993 November 11 flare and 1992 October 5 flare, 
respectively. They found the value of 0.003 - 0.2 in preflare and impulsive phase, and 0.0002 
- 0.0003 in gradual phase. Isobe et al. (2002b) analyzed the decay phase of an LDE flare on 
1997 May 12 and estimated the reconnection rate at 0.001-0.01 using the energy release rate 
and the reconnected magnetic flux. Using the same method as the previous one, Isobe et al. 
(2002a) analyzed 7 flares and 2 giant arcades and found that the reconnection rate is 0.001 
-0.1. Moreover, analyzing the impulsive phase of three two-ribbon flares, Isobe et al. (2005) 
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obtained the value of 0.015 - 0.071. 

Recently, direct measurements of the reconnection inflow were reported. Yokoyama 
et al. (2001) found an evidence of reconnection inflow and derived the reconnection rate as 
0.001 - 0.03. By the same method, Narukage & Shibata (2006) reported 6 inflow events and 
found the reconnection rate is 0.001 - 0.07. Lin et al. (2005) studied an eruptive limb event 
occurred on 2003 November 18 and reported the reconnection rate is in the range from 0.01 
to 0.23. Note also that they found the reconnection inflow using the Doppler shift. 

In this paper, we report a statistical study of flares that occurred in the year of 2000 
observed with Yohkoh/SXT. We determine several physical quantities of each flare and try 
to estimate the reconnection rate in solar flares. Our aim is to find statistical tendency of 
the physical parameters. In §2, we describe the data and the method of analysis. In §3, we 
show results. Discussion and conclusion are given in §4 and §5, respectively. 



DATA ANALYSIS 



2.1. Inflow Velocity 



The amount of energy released during a flare Ea arc can be explained by the magnetic 
energy stored in the solar atmosphere; 
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where L is the characteristic size of the flare and B COT is the characteristic magnetic flux 
density in the corona. Since the released magnetic energy balances the energy flowing into 
the reconnection region, we can describe the energy release rate as 
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where V in is the inflow velocity of the plasma. Therefore, the time required for the energy 
inflow to supply the flare energy is estimated as 
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and this should be the timescale of the flare. Using this timescale, we can estimate the inflow 
velocity V in as 
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To evaluate the reconnection rate in non-dimensional form Ma = V- ni /Vx, we must estimate 
the Alfven velocity in the inflow region Va = B COT / ^Airp. Hence, if we measure the coronal 
density p, the spatial scale of the flare L, the magnetic flux density in the corona -B C or? and 
the timescale of flares Tfl arc , we can calculate inflow velocity V in , Alfven velocity Va, and 
reconnection rate Ma. 

2.2. How to Determine the Physical Parameters of the Solar Flare 

In this subsection we describe how to determine the parameters required for derivation 
of the inflow velocity V; n and the Alfven velocity Va- To estimate the inflow velocity V; n by 
using equation (4), we need the spatial and temporal scales of each flare: L and Tfl are . First, 
the flare size L is determined by soft X-ray images taken by the Soft X-ray Telescope (SXT; 
Tsuneta et al. 1991) onboard the Yohkoh satellite (Ogawara et al. 1991) with the beryllium 
(Be) filter. We define L as the square root of the area where the intensity is more than 10% of 
the maximum value in each frame and it is measured at the time when GOES X-ray intensity 
attains its maximum value in each event. The error of L is considered as the difference in 
definition of the size; the minimum size is defined by the area where the intensity is more 
than 30% of the maximum value in each frame and the maximum size is defined by the area 
where the intensity is more than 3% of the maximum value. The timescale of a flare Tfl are is 
defined as the period from the onset to the peak of GOES X-ray flux 1 . Since we use one-min 
average of GOES X-ray flux data, the error of this timescale is considered as one minute. 
In equations (1) and (2), the spatial scale of the flare L should be the length of the current 
sheet. Although we cannot measure it or estimate the ratio of the current sheet to the global 
size accurately, we can consider the length of the current sheet is the same order as the size 
of system since the flare size L is measured in the impulsive phase. 

To estimate the Alfven velocity in the inflow region Va, we need coronal density p and 
coronal magnetic flux density B cor . Owing to difficulty in measuring the coronal density p, 
we assume it 1.67 x 10~ 15 g cm -3 that is the typical value in the active region. The coronal 
magnetic flux density B cor is also difficult to measure directly. Therefore, we try to evaluate it 
in two ways. In the first method (method 1), we use the photospheric magnetograms from the 
Michelson Doppler Imager (MDI; Scherrer et al. 1995) onboard the Solar and Heliospheric 
Observatory (SOHO; Domingo et al. 1995). We define the representative magnetic flux 
density for each flare as follows: First, using soft X-ray images taken by SXT with the Be 
filter, we determine the area in which the intensity is larger than a threshold value. Second, 
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we measure the magnetic flux in the area and obtain the average of unsigned magnetic flux 
density. Changing this threshold value from 10% to 50% of the maximum soft X-ray intensity 
of each image, magnetic flux densities are obtained as a function of the threshold value of 
the soft X-ray intensity Finally, we define representative magnetic flux density as the 
midpoint of the maximum and the minimum of the obtained flux density in this range and 
consider the extent from the maximum to the minimum as the error range of _B p h. 

In Figure 1 we show the result of the GOES M3.9 class flare on January 18 as an 
example. The left panel of Figure 1 illustrates the photospheric magnetogram of the flare 
region. White and black indicate positive and negative polarities in the magnetogram. 
The overlaid contours correspond to 10, 30, and 50% of the maximum soft X-ray intensity. 
The right panel of Figure 1 shows photospheric magnetic flux density as a function of the 
threshold value of the soft X-ray intensity. The dotted lines indicate the maximum and 
minimum flux density and the solid line corresponds to the midpoint value, which is defined 
as the representative photospheric flux density of this event. 

In order to obtain the coronal magnetic flux density, we assume 

B COT = a£5 ph (5) 

and 

a B = const. = 0.3 (6) 
among all flares. The validity of this assumption will be discussed in §4.1. 

The second method (method 2) to obtain the coronal flux density is as follows. If all the 
magnetic energy stored in the corona before the reconnection process turns into the thermal 
energy in the reconnected loop, we obtain 

B 2 

Ploop ~ g^> ( 7 ) 

where pi oop is the plasma pressure in the flare loop and B eq is the magnetic flux density in 
the corona. To distinguish this magnetic flux density from that derived by method 1, we use 
the subscript "eq" . The value of j9i oop is estimated from the color temperature T and the 
emission measure EM calculated from the GOES X-ray data (Garcia 1994), that is, 

Pioop = ni oop k B T, (8) 

where 

I EM 

"-loop ~ 2 y-^T' ( 9 ) 

ni 00 p is the number density in the flare loop, and k-Q is the Boltzmann constant. Since 
the GOES X-ray flux data have no spatial resolution, we use the flare size measured in 
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the SXT image as the size L in equation (9). Hence, we can estimate B cq and assume 
B cor ~ B cq . The error of this flux density B eq results from the error in deriving temperature 
and emission measure from the GOES X-ray data. On the basis of the description in Garcia 
(1994), estimated uncertainty in the derived temperature depends on the temperature and is 
calculated at 13% around 10 MK. We calculate this uncertainty as a function of temperature 
and then obtain the propagated error of B cq . 



2.3. Event Selection 

As we describe in the previous subsection, we use SXT images, MDI magnetograms, and 
GOES X-ray flux data. We survey flares that occurred in 2000 under the following criteria: 
(1) GOES class of the flare is above C6.0, (2) the flare site is included in the field of view of 
the partial-frame SXT image with the Be filter taken at the peak of the GOES soft X-ray 
flux, and (3) the distance from the disk center to the flare region is within 800 arcsecs to 
avoid the line-of-sight effect in estimation of the magnetic flux density. From the flare list 
obtained by GOES' data, 482 flares meet the above criterion (1) and 77 of those flares satisfy 
all the criteria. 



3. RESULT 

Using the method described in the previous section, we have analyzed 77 flares occurred 
in the year of 2000. First, we examine dependence of physical parameters of these flares on 
GOES class. In Figure 2 we plot the parameters of flares as a function of GOES class. 
Temporal scale Tfl arc , spatial scale L, photospheric field £> p h obtained by method 1, and 
coronal field B eq obtained by method 2 are shown in Figure 2a, 2b, 2c, and 2d, respectively. 
We find from these plots: (1) There is weak correlation between the timescale Tfl arc and the 
GOES peak flux (Panel a). This tendency is consistent with Figure 3 of Veronig et al. (2002) 
and Figure 2 of Kay et al. (2003). (2) The characteristic size of flares L shows a larger scatter 
when the GOES' peak flux is smaller (Panel b). Garcia (1998) reported the relation between 
the GOES peak flux and the loop length of flares in their Figure 5. Comparing our results 
with theirs, the value of the size is in the same range, although clear correlation between 
these parameters cannot be found out in their results. (3) There is a threshold value in the 
magnetic flux density (£> p h and B cq ) that increases as the GOES peak flux increases (Panels 
c and d). For example, the X-class flares occurred only when B p ^ was larger than 100 
Gauss. 



-7- 



Second, we examine the relationship between the physical parameters of flares. Figure 
3 shows the spatial scale L plotted against the timescale Tfl are and tells that the spatial scale 
L tends to be larger with increasing timescale. This result is consistent with the result 
shown in Figure lc of Garcia (1998), although the tendency is not so clear in their figure 
due to a small number of their analyzed events. The lines in Figure 3 indicate the inflow 
velocity V in given by equation (4). The inflow velocities are found to be distributed from a 
few km s _1 to several tens km s _1 . These inflow velocities are comparable to those in the 
previous studies. Detailed discussion is given in §4.2. The left panel of Figure 4 shows the 
photospheric flux density 5 p h obtained by method 1 plotted against the timescale Tfl are , and 
in the right panel of Figure 4 the coronal flux density B cq obtained by method 2 is plotted 
against the timescale Tfl are . Assuming the ratio of coronal flux density to photospheric flux 
density «g = B^/B^ is 0.3 and the coronal density p is 1.67 x 10 -15 g cm~ 3 , the values 
of magnetic flux density in the left panel of Figure 4 correspond to 10 3 - 10 4 km s _1 of the 
Alfven velocity in the inflow region Va- The values of B eq in the right panel of Figure 4 also 
correspond to the same range of V\. 

Using these results of Vi n and Va, we examine the reconnection rate Ma- In the left 
panel of Figure 5, the reconnection rate obtained by method 1 is plotted against the GOES 
class. The values of reconnection rate are distributed around 10~ 3 - 10~ 2 . Note again 
that we assume = i? cor /-Bph = 0.3 in this plot. By changing as, the data points will 
shift proportionally to the value. The right panel of Figure 5 shows the reconnection rate 
calculated by using the magnetic flux density B eq obtained by method 2 against the GOES 
class of each flare. The reconnection rates are also distributed around 10~ 3 - 10~ 2 . Both 
panels of Figure 5 show that the reconnection rate tends to diminish as the GOES class of 
the flare increases. 

As mentioned in §2.2, we estimate the coronal magnetic flux density in two ways. In 
Figures 2c and 2d, although the flux density obtained by method 1 (Figure 2c) appears more 
broadly scattered than that obtained by method 2 (Figure 2d), the values of flux density 
show the similar distribution. Moreover, the magnetic flux density obtained by these two 
methods has a similar distribution as a function of Tfj are (Figure 4). Two panels of Figure 
5 show that the values of the reconnection rate obtained by methods 1 and 2 are the same 
range and display the same dependence on the GOES class. Hence, we conclude that the 
results by these two methods are consistent with each other. 
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4. 



DISCUSSION 



4.1. Ratio of Coronal Flux Density to Photospheric Flux Density 



We must examine whether the assumption that the ratio of coronal flux density to 
photospheric flux density a B = B cor /B ph equals 0.3 is appropriate one. In Figure 6, we 
compare the photospheric flux density B p \ l obtained by method 1 with the equipartitional 
flux density B eq obtained by method 2. The lines in the figure correspond to 0.1, 0.3, and 1 
of Bcq/Bph. This figure shows the ratios are distributed around 0.3, and if we consider B eq 
as the coronal flux density, our assumption a B = 0.3 is not inappropriate. 

Dere (1996) estimated coronal flux density from data of photospheric flux density de- 
scribed in the scientific literature. He calculated typical magnetic flux density as a function 
of height as photospheric magnetogram forms the base. The formula for his analyzed active 
region observed on 1974 January 16 is 



where h is height above the photosphere. He applied this formula to calculate the coronal 
flux density of all his analyzed events in active regions. Using the typical spatial scale of 
flare loops in our analysis L ~ 10 9 cm as the height h, we obtain the coronal field of 38 
G by equation (10); we can also calculate the photospheric flux density of his model at 95 
G by taking h = in equation (10). Therefore the ratio of the coronal flux density to the 
photospheric flux density is roughly estimated at 0.4. This formula probably underestimates 
the magnetic flux density, because the photospheric flux density analyzed in our paper is in 
the range from 100 G to 700 G. Hence, we also did a calculation of the coronal potential field 
for some particular events for more investigations. The left panel of Figure 7 shows average 
potential magnetic field strength of X2.3 flare occurred on 2000 November 24. This flux 
density is extrapolated by using the software package MAGPACK2 (Sakurai 1982) and its 
strength is determined and averaged in the central portion of the flare. In this calculation, 
the coronal flux density at the height of 2.56 x 10 9 cm that corresponds to the flare size of 
this event is 83 G. As shown in the left panel of Figure 7, the photospheric field strength is 
335 G at h — 0. Note, however, that as in equation (5) should be calculated by the ratio 
of the coronal flux density to the line-of-sight strength of the photospheric field. Therefore 
we calculate the line-of-sight component of photospheric flux density using calculation of the 
potential field; the spatial average of the line-of-sight component at h = in the flare site 
is B ph = 151 G, which is shown as an asterisk in the left panel of Figure 7. Accordingly, 
we obtain the ratio of coronal flux density to the line-of-sight component of photospheric 
flux density as is nearly 0.6. In the right panel of Figure 7, we show the average potential 
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magnetic field strength of M3.7 flare occurred on 2000 July 14 in the same way. Concerning 
this event, the photospheric flux density is 214 G at h — 0, the coronal flux density at the 
height of the flare loop is 22 G and the average line-of-sight component of photospheric field 
is B ph = 119 G; is roughly 0.2. We calculate some other events and obtain the ratio «b 
in the range of 0.2 - 0.6. 

Considering these estimates of the ratio, the assumption a B e B cor /B ph ~ 0.3 is proba- 
bly plausible. As a result, on the basis of the result described in §3, we find the reconnection 
rate Ma is 10~ 3 for our analyzed events. 

4.2. Comparison with the Previous Study 

In this section, we compare our results with the previous studies to evaluate our results. 
Using inflow pattern in extreme ultraviolet images, the inflow velocity V; n was measured 
around 5 km s _1 by Yokoyama et al. (2001), while Narukage & Shibata (2006) obtained 2.6 
- 38 km s _1 in six events. Including other indirect measurements, the inflow velocity is in 
the same range as our results (Figure 3). Moreover, we analyze one of Narukage & Shibata's 
events on 1999 May 27 using our method described in §2.2. We estimate the inflow velocity 
at 5.2 km s" 1 , while they reported the inflow velocity of this event is 12.5 - 14.9 km s" 1 or 
31.3 - 37.3 km s _1 . The former is obtained by the method based on Yokoyama et al. (2001), 
while the latter is obtained by the estimation based on Chen et al. (2004). In either case, 
the velocity obtained by our method is less than that obtained by their method. 

Next, we compare our results with those of Isobe et al. (2005), since some of our analyzed 
events are also examined by them. In our statistical study we use relatively simple methods 
to study a large number of events, while they dealt with three events in detail in their case 
study. Table 1 summarizes parameters of the flares obtained in this work and in Isobe et al. 
(2005), mainly referring to their Table 1. We calculate some parameters based on this table 
as follows. Isobe et al. (2005) measured the width L x , the length L y , and the height L z of 
the flare arcade separately. To compare with our measurement of L, we calculate the flare 
size L is defined as the cubic root of the product of L x ,L y , and L z given in the paper. The 
averaged energy release rate E^^/r is defined as the thermal energy in the flare loop divided 
by the temporal scale, 

3 

Efi^e/r— -ri\ oop k B TL /r flare . (11) 

This value is calculated from the peak value of E t ^ in Isobe et al.'s Figures 6, 8, and 10 to- 
gether with our GO.SS'timescale Tfl are . Moreover, we define the energy release rate \dE mag /dt\ 
as equation (2) for our study and as H in Isobe et al. (2005) for their study, both of which 
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correspond to the Poynting flux into the diffusion region. Furthermore it should be noted 
here that the GOES C-class event in Table 1 is not one of our analyzed events in this statis- 
tical study due to lack of the non-saturated SXT image with the Be filter at the flare peak. 
Therefore using the last non-saturated image with the Be filter taken at three minutes before 
the peak, we calculate the parameters in this table for only comparison. 

Comparing our result with Isobe et al. (2005), we find that in our study the inflow 
velocity V[ n is systematically smaller. Since we calculate V- Ul as L/(4rfl are ), V- in obtained in 
our study is an average value in the impulsive phase of a flare. On the other hand, Isobe 
et al. (2005) estimated Vj n using the peak value of the energy release rate in the impulsive 
phase of the flares. Therefore our small inflow velocity Vj n perhaps reflects nonsteadiness 
of the reconnection rate in the reconnection process. Moreover, owing to our larger coronal 
magnetic flux density B COT , the coronal Alfven velocity V\ of our result is systematically 
larger than those of their result. Consequently, this larger Va and smaller V in make the 
reconnection rate Ma smaller than a factor of 2 - 10. Nevertheless, considering energy 
budget of each flare, we note that our result is not inconsistent; the energy release rate 
\dE mag /dt\ is comparable with the averaged energy release rate -Efl are /r during the impulsive 
phase of the flare. 



4.3. Reconnection Rate and Reconnection Model 

Here we discuss the possible implication to the models from the reconnection rate ob- 
tained in this work. As mentioned in §1, different models show not only different values of 
the reconnection rate but also different dependence of the reconnection rate on the magnetic 
Reynolds number R m = V^L/i]. Therefore we examine this dependence. The reconnection 
rate Ma is plotted against the magnetic Reynolds number R m in Figure 8. The reconnec- 
tion rate and the Reynolds number obtained by method 1 is shown in the left panel of the 
figure and those obtained by method 2 is shown in the right panel. In the left panel, the 
dash-dotted line indicates the least-squares fitting of the data, which yields M A oc -R~ ' 8 , 
while in the right panel there is no clear dependence since the magnetic Reynolds number 
is distributed within one order of magnitude. To see the figure clearly, the error bars of the 
magnetic Reynolds number are omitted in the right panel, though the error is one order of 
magnitude at most. The dotted lines in both panels indicate the dependence of the recon- 
nection rate on the magnetic Reynolds number that is predicted by the Petschek model and 
the dashed lines indicate the Sweet-Parker-type dependence Ma oc -R m ' 5 - In calculation of 
the magnetic diffusivity i] in R m , we assume the resistivity is attributed to the Coulomb 
collision (Spitzer 1956) and the plasma temperature is 10 7 K. 
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In the process of estimating the reconnection rate Ma and the magnetic Reynolds 
number R m shown in Figure 8, we have following assumptions: (1) the coronal density 
p is 1.67 x 10~ 15 g cm -3 , (2) the ratio of coronal flux density to photospheric flux density 
«b = B cor /B p h equals 0.3, and (3) the plasma temperature is 10 7 K. If the coronal density 
p and/or the ratio ats deviate from the assumed value, the values of the Alfven velocity and 
the magnetic Reynolds number are affected. However these values will only shift by some 
factors on the whole. Therefore the distribution shown in Figure 8, i.e., the power index 
of the dependence of the reconnection rate M A on the magnetic Reynolds number i? m will 
not be greatly affected by assumptions (1) and (2). Moreover, we can consider the factor 
to shift to be less than one order of magnitude, because the coronal density p is thought to 
be within at most one order of magnitude from the typical value and the ratio «g is esti- 
mated in the range of 0.2 - 0.6 as mentioned in §4.1. Furthermore, in analyzing by method 
2, we calculate the plasma temperature from GOES X-ray data and find that the observed 
temperature of plasma is distributed within less than one order of magnitude around 10 7 K. 
Hence assumption (3) is thought to have little effect on our results. 

From the comparison of ours with the previous results mentioned in §4.2 (e.g., Isobe 
et al. 2005; Narukage & Shibata 2006), the error of the obtained reconnection rate shown 
in our figures is probably a little underestimated. Besides the assumptions mentioned in 
the previous paragraph, other uncertainties, such as the uncertainties in the measurements 
of the size L, the timescale Tfl are , the temperature T and the emission measure EM, will 
affect our result. Even after considering all these effects, the value of M A in Figure 8 is 
found within one order of magnitude from the theoretical value of the Petschek model, 
while dependence of the reconnection rate on the magnetic Reynolds number tends to be 
stronger than that in the Petschek model. This result suggests two possibilities: One is 
the occurrence of the Petschek-type reconnection in the flares and the other is the Sweet- 
Parker model with the magnetic diffusivity rj enlarged by some diffusion processes, such 
as the MHD turbulence. Owing to the uncertainties, we cannot immediately rule out the 
possibility of either the Petschek model or the generalized Sweet-Parker model from our 
result. We expect to discriminate between these two interpretations by the observations of 
Solar-B. If the Sweet-Parker reconnection with this enhanced resistivity occurs in a flare, the 
thickness of the diffusion region is estimated to be 100 km, when the length of the diffusion 
region is approximately the size of flare L ~ 10 9 cm and the reconnection rate is Ma ~ 0.01. 
Because the size of diffusion region is less than the limit of resolution, the structure of 
turbulence is expected to be observed as a single-pixel line-broadening with Solar-B. On the 
other hand, if the Petschek-type reconnection occurs, the slow shock structure is expected 
to be observed. 
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5. CONCLUSION 

We analyze the data of flares occurred in the year of 2000 and measure the temporal 
scale r, the size L, and the magnetic field B of each flare. First, we examine the dependence 
of these parameters on the GOES class and find that the sizes of flares tend to be distributed 
more broadly as the GOES class becomes weaker, and there is a threshold value in both pho- 
tospheric and coronal magnetic flux density that increases as the GOES peak flux increases. 
Second, we examine the relationship between those parameters. There is weak correlation 
between temporal and spatial scales of flares, while there is little correlation between tempo- 
ral scale and magnetic flux density. Finally, we try to estimate the reconnection rate in each 
flare. The values are distributed in the range from 10~ 3 to 10" 2 and the reconnection rate 
decreases as the GOES class increases. The value of the reconnection rate obtained in our 
study is within one order of magnitude from the predicted maximum value of the Petschek 
model, although the dependence of the reconnection rate on the magnetic Reynolds number 
tends to be stronger than that in the Petschek model. 
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erated by the National Oceanic and Atmospheric Administration's (NOAA's) National En- 
vironmental Satellite, Data, and Information Service (NESDIS). SOHO is a project operated 
by the European Space Agency and the US National Aeronautics and Space Administration. 
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Table 1. Parameters of the Flares 



Parameter 


2000 Nov 24 15:13 X2.3 flare 


2000 Jul 14 13:52 M3.7 flare 


2000 Nov 16 0:40 C8.9 flare 




This Work 


Isobe et al. (2005) 


This Work 


Isobe et al. (2005) 


This Work 


Isobe et al. (2005) 




Method 1 Method 2 




Method 1 Method 2 




Method 1 Method 2 




T (sec) 


1320 




480 




1200 




Bph(G) 


205 ± 17 


449 


106 ± 11 


117 


30 ±4 


106 


Bcor(G) 


62 116 


41 


32 60 


44 


9 32 


11 


L(10 9 cm) 


2.56 


2.8 


2.94 


2.3 


4.12 


4.0 


T(10 6 K) 


14 


10 


8.8 


8.0 


7.2 


8.0 


y in (cm s- 1 ) 


4.8 x 10 5 


1.3 x 10 7 


1.5 x 10 6 


3.2 x 10 6 


8.6 x 10 5 


6.7 x 10 6 


V A (cm s- 1 ) 


4.25 x 10 s 8.0 x 10 8 


2.8 xlO 8 


2.2 x 10 8 4.2 x 10 8 


2.1 x 10 8 


6.2 x 10 7 2.3 x 10 8 


9.4 x 10 7 


M A 


1.1 x 10~ 3 6.0 x 10~ 4 


4.7 x 10~ 2 


7.0 x 10~ 3 3.7 x 10" 3 


1.5 x 10~ 2 


1.4 x 10~ 2 3.9 x 10~ 3 


7.1 x 10~ 2 


EflareMerg s" 1 ) 


1 x 10 28 


4 x 10 27 


1 x 10 28 


3 x 10 27 


4 x 10 27 


8 x 10 26 


|<2£ mag M|(erg s" 1 ) 


2 x 10 27 7 x 10 27 


2.7 x 10 2Sa 


2 x 10 27 8 x 10 27 


6.4 x 10 27a 


2 x 10 26 2 x 10 27 


1.9 x 10 27a 



a Variablc H in Isobe et al. (2005) 

I— > 
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X (arcsecs) 



Fig. 1. — Left panel: Soft X-ray contours overlaid on the SOHO/MDl magnetogram of the 
M3.9 flare on 2000 January 18. The contours show Yohkoh/SXT image; 10, 30, and 50% of 
the maximum soft X-ray intensity. White and black indicate positive and negative polarities 
in the magnetogram. Right panel: Photospheric flux density plotted against threshold value 
(see text). The dotted lines indicate the maximum and minimum of the flux density and the 
solid line corresponds to the midpoint of the maximum and the minimum. We define this 
midpoint as the representative photospheric flux density of this event and the error range as 
the extent from the minimum to the maximum. 
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Fig. 2. — Physical parameters of each flare plotted against the GOES class, (a) Time scale 
Tflarc- (b) Size L. (c) Photospheric magnetic flux density B p \ l obtained by method 1. (d) 
Coronal flux density B eq obtained by method 2. 
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Fig. 3. — Spatial scale of each event plotted against timescale. Crosses, triangles, and squares 
indicate X-class, M-class, and C-class, respectively. The solid line, dotted line, dashed line, 
and dash-dotted line indicate 3, 10, 30, and 100 km s _1 of the inflow velocity V in , respectively. 
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Fig. 4. — Left panel: Photospheric flux density of each event B ph , obtained by method 1, 
plotted against timescale of the flares Tfl are . The dashed lines indicate the corresponding 
coronal Alfven velocity calculated on the assumption that the coronal density p = 1.67 x 
10~ 15 g cm -3 and a# = B COI / B^ = 0.3. Right panel: Coronal flux density B eq , obtained by 
method 2, plotted against the timescale r flare . The dashed lines indicate the corresponding 
coronal Alfven velocity. Meaning of the symbols in these figures is the same as Figure 3. 
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Fig. 5. — Reconnection rate M A plotted against the GOES class of each flare. The left panel 
shows the reconnection rate obtained by method 1. We assume «b = -B cor /-Bph = 0.3 when 
we calculate Ma- The right panel shows the reconnection rate obtained by method 2. 
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Fig. 6. — Coronal magnetic flux density B eq plotted against photospheric flux density B ph . 
The solid line, the dotted line, and the dashed line correspond to 0.1, 0.3, and 1 of B eq / B ph , 
respectively. Meaning of the symbols is the same as Figure 3. 
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Fig. 7. — Spatial average of potential-field flux density plotted as a function of height. The 
left panel shows an example of 2000 November 24 X2.3 flare and the other example of 2000 
July 14 M3.7 flare is shown in the right panel. The dashed line indicates the size of each 
flare and the asterisk represents the line-of-sight component of the photospheric flux density. 
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Fig. 8. — Reconnection rate Ma plotted against the magnetic Reynolds number R m . The 
left panel shows the reconnection rate obtained by method 1, and the right panel shows 
the reconnection rate obtained by method 2. In both panels, the dotted lines indicate 
the dependence of reconnection rate on the magnetic Reynolds number that is predicted 
by the Petschek model, and the dashed lines indicate the Sweet-Parker-type dependence 
M\ (x -R^ ' 5 - The dash-dotted line in the left panel indicates the least-squares fitting of the 
data, which yields Ma oc -R^ ' 8 - Meaning of the symbols is the same as Figure 3. 



